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ABSTRACT 



The Haltlner theory of saturated steady state 
convective currents which includes the transfer of heat 
and momentum by lateral diffusion as well as the 
entrainment of environmental air is discussed. A 
modified model is presented which includes, as a 
separate term, the drag of the liquid water carried 
along by the ascending current. The system of equations 
is integrated numerically for several cases previously 
considered by Haltlner allov/ing entrainment to assume 
negative values (detrainment) and the results are 
compared. The results are in good agreement, and show 
that the net effect of the opposing actions of liquid 
water drag and detrainment is a reduction in the 
computed vertical velocities. This reduction appears 
to be a correction in the right direction on the 
basis of a comparison to some observed data. 

The writer wishes to express his appreciation for 
the assistance and encouragement given him by 
Professor G.J. Haltlner of the U.S. Naval Postgraduate 
school in this investigation. 
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Introduction 



The concept of simple "adiabatic parcel ascent" as 
a first approximation to convective phenomena has long 
been in use; it has, however, a number of shortcomings. 
As a result of evidence (Stommel,[g] ) that a rising- 
convective current acts like a jet, dragging in surround- 
ing air as it progresses upward, Austin and Fleisher 
suggested that this process, called "dynamic entrain- 
ment", is necessary to satisfy mass continuity. t.l] 
HaltinerC^} summarized the work of Houghton and Cramer, 
Bunker, [2I and MalkusjfS'J and suggested a theory of 
saturated convective currents which includes the effects 
of lateral diffusion of heat and momentum, and 
extended computations to elevations typical of large 
cumulus clouds. 

Haltiner [Hj derived a system of equations which 
describe a saturated convective current. The derivation 
of these equations is essentially the same as the 
derivation of the modified system of equations presented 
below. The Haltiner model assumes that during 
detralnment the cross-sectional area of the rising 
column remains constant. However, when mass continuity 
considerations calls for the ejection of mass 
(detralnment), a second system of equations is applied 
based on the assumptions that the rate of entrainment 
remains identically zero and mass continuity is 
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satisfied by an increasing cross-sectional area in 
the rising air. In this paper the detrainment phase 
will be treated by retaining the assumption that the 
cross-sectional area of the rising air remains constant, 
thus permitting the entrainment to assume negative 
values o 

The effect of the drag of the liquid water carried 
along in the rising air was not explicitly considered 
in the Haltlner model, although it can be qualitatively 
Included by adjusting the rate of diffusion of 
momentum. Here the drag effect of the liquid water will 
be specifically included as a separate term. 

The following notation will be used throughout 
this papers 
T temperature 

z height 

q specific humidity 

Cp specific heat of air at constant pressure 

L latent heat of vaporization 

g gravity 

w velocity 

k diffusion coefficient 

gas constant for dry air 
R,^ gas constant for water vapor 

N In M 

M Mass rate of flow per unit time 
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j( "specific humidity" of liquid water 
ff cross-sectional area of rising column 



density 



p pressure 
lapse rate 



subscripts : 
e environment 
V virtual 
o initial value 
2o The system of equations 

a» The continuity equation. 

Consider a vertical current of air in which the velocity, 
w, is horlzontall”’ uniform. If the field of vertical 
velocity is not uniform w may be considered to be the 
mean value in the horizontal at any given level. The 
mass of air per unit time flowing past any given level, z , 



where M is the mass,/> the density and (T the cross- 
sectional area of the vertical current. Under steady 
state conditions the mass at an adjacent level Z +dzls 
the sum of the mass rate of flow at level 2 , and air 
entrainment from 2 to z + dz • Thus, 

ol i'^ cl^ cl w ci (T 



is given by 



/v\ vv a 



-V 
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( 2 ) 



Defining a parameter N ■= InM, we may write 

cl N ^ cj^ ^ < dw J_ J/;- 
d X ol z w ci i (T cl z 



Differentiation of the equation of state p - 
yields 

/ox 1 = -Lfif - — — ^ 

/"dz pdz 

We next make the assumption that at any level the 

pressure in the rising air is the same as that of the 

environment, then _L =■ _L c)f^ 

P dz (3 

From the hydrostatic equation ~ 9^^ 

follows that 

(4) 1 ie = = - J— 

P 1’ Ra'^^e 

We may now eliminate the density term from (2) by 
substitution from (3) and (4) to yield the continuity 
equation for the rising air 

( 5) dN _ j_ civv ^ _ J_ dX 

dz w dz (T dz X 

bo The Equation of Motion. 

The changing mass of the rising air due to entrainment 
or detrainment requires the use of an equation of 
motion for a variable mass. If a mass M is moving at 
velocity w at time t combines with a second mass dM 
moving at velocity w^ during a time interval dt 
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and w dw represents the velocity of the combined masses 
at time t + dt, then the momentums at time t and t t dt 
are related by the equation 

( M +c1M)(vv +d^) Mw" + t M cjt 

Here E represents the external forces per unit mass. 

This relationship yields the equation of motion 



(6) ^ -t- M 

clt M d t 









J. 

I'l clt 



The external forces acting on the system include the 
pressure and gravity forces, which may be expressed as 
g { X" frictional forces P, 

the latter including the drag of the liquid water carried 
by the ascending column of air. The equation of motion 
for the liquid water rising with the air may be written 



( 7 ) J ^ _ J ^ -V- D 

d t 

where is the mass of liquid water per unit mass of 
air which we will term the "specific humidity" of the 
water and D represents the drag of the air acting on 
the liquid water of massj^ . The negative of D is the 
drag of the liquid water on a unit mass of air. From 

eq (7), we have 

-0 = - 

substituting in eq (6) we obtain 

ire “ 9 ( j - j id H- F 

^ * V<2 M <3 "t 

simplifying we have 

+ p 

® i v-e, t 
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assuming steady state conditions we obtain 





Co Lateral diffusion of heat and momentum 
In addition to the exchange due to entrainment, 
turbulent exchange of heat and momentum must also be 
consideredo Priestley C73 ( 1956) investigated the 
mechanisms of diffusive heat and momentum exchange in 
a plume of hot gases » Haltlner [^] suggested that, 
except for the scale of the phenomena, the diffusive 
exchange processes in a convective current are similar 
to the processes defined by Priestley. The effects of 
lateral diffusion were included in the Haltlner model 
by the extention of the Priestley equations to convective 
currents. The heat loss per unit time is taken to be 



It is also necessary to saturate the environmental air 
entrained into the updraft and a term of the form 



is required, where L is the latent heat of vaporization. 
The loss of momentum per unit time per unit mass is 



( 9 ) 




or approximately 



(10) -A, c,, (t-t;) 



(11) L (a- 4c) 
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similar to that of the loss of heat and is assumed to 
be proportional to the difference between the velocity 
of the convective current and the environment, with the 
latter taken to be zero, with the result — where 
has the dimension second”,* W ith this form for the 
friction term, eq (8) becomes 



(12) 



cl 



>A/ 



o s 



(i^i) 



9^ 

w 



( X~ X a. _x] - (. W - Wg ) 4S. - V?, 

^ X. ^ a z 



The dimensions of imply that diffusive losses are 
independent of updraft velocity. Haltiner suggests a 
second approach to the problem of diffusion which makes 
the loss of heat and momentum dependent upon updraft 
velocity. The first approach, however, will be employed 
in this paper. 



d. Thermodynamic equation. 

In describing the thermodynamic process undergone 
by the rising air a number o'^ factors must be included. 
The following must be taken into account in accordance 
with the First Law of Thermodynamics; the heat released 
by condensation, the heat required to saturate the 
entrained environmental air and to raise its temperature 
to that of the rising column, the heat lost to the 
environment through diffusion, the internal energy 
changes and the work done during expansion. The 
resulting equation embodying these factors is as 
follows; 
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(]_3) —MLidc[ — L(c[-<^e)^l 

- /^Jk, -M_c^--4j]cit - 

/V) [c^,clT- 

The terms on the left, in the order of their 
occurrence, represent respectively, (a) the heat 
released In condensation; (b) the heat required to 
saturate the entrained air, (c) the heat required to 
bring the entrained air to the temperature of the rising 
column, (d) the heat lost through diffusion in the 
interval dt required for the air to rise a distance 
dZo The right side of the equation, represents the 
combined internal energy change and work by the gas in 
expanding from pressure p to pressure p-vdp. By 
substitution from eq (4), eq (13) becomes, with 
appropriate manipulation. 



- '-‘ii - le)-vL(ci-c^)]cil^ 

az dz 






= c.,OT 



ciZ Xe 
e. Moisture and liquid water. 

It has been assumed that the rising air is sat- 
urated; the specific humidity of the rising air is then 
a function of temperature and pressure. As a good 
approximation we have the expression = .62.l‘^p where e 
is the saturation vapor pressure. By differentiation 
this becomes 
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(15)j_ clq _ j_ _ 1 iP 

q dz ' ^ cU p clz 

The last terra in eq (15) ^nay be replaced with 

from eq (9). and from Clapeyron's equation, — - — ^ , 

j (± i ^ ’ 

y 

where is the gas constant for water vapor, an 
alternate expression for the second terra is obtained, 
yielding 



( 16 ) 



I del _ Jz sL! 



A 



R cl '2 



RT-.b R.T, 



<i V e 



The liquid water content of the cloud is the amount of 
moisture condensed less the amount required to saturate 
the environmental air entrained and that exchanged 
through diffusion. Defining J* as the "specific humidity" 
of the liquid water we may write 

(1?) d(nt) =- hclc^ -(4-%) clK\ - 



Expanding eq (1?) yields 



(18),di+ ^ dN 

dz clz <dz 



+ Al =0 

VV 



Equation (17), however, is not valid during detrainment 
and the change in liquid water content must be piece- 
wise defined. For the detrainment phase we have 



(19) cl (Mi) = - hd<a + Hd M - tnk, (c[-qe) at 

which upon expansion leads to 
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(20) _ _ ii _ Jii 

cl z <ci K 

Two additional relationships between temperature and 
virtual temperature will be required. One of these is 
the well known relationship 



(21) ~T — ( I i- . 6 nr 



the second relationship is obtained from eq ("I) by 
differentiating with respect to z ; this yields 
^ = (it ^ 

^ <J ^ 

Substituting from eq (l6 ) this may be written 

— = V R„T Xlz Pf "T 



R,T 

<i ve 



J Z 

By using some mean values this may be approximated by 



( 22 ) d~R _ /l + 

d ^ " i /Jz 

f. The environment. 

The environmental parameters , and 

appear in one or more of the equations derived above; 
it is therefore necessary to define the environment in 
some way. We have wide latitude in choosing the struc- 
ture of the environment; it may, for instance, conform 
to an actual sounding. A simple model was used in this 
investigation. The lapse rate and the relative 
humidity were assumed constant, and w^ was assumed to 
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be zero. Under these assumptions may be expressed 



as 




The environmental specific humidity may be determined 
as follows: From eq (l6) the saturation specific 
humidity is given by 



and since the relative humidity is constant we 

may write 



If the lapse' rate , the initial temperature"!^^ , and 
the relative humidity are chosen, then the values ofT^,, 

, and may be obtained at any level from equations 
(21), ( 23 ), and (2^). 

3. The complete model. 

Having defined the environment in sub sc-;:, vic'n r 
equations (5)> (12), (14), (l6), (21) and (18) or (20) 
constitute a system of six equations in 7 unknowns, w, 

T ,'1^^ >q > i > H, and (T . We must now establish a 
relationship defining <f . The simplest solution, and 
the one which will be employed here, is to assume that 
(T is Identically zero. This assumption is quite 
likely an oversimplification, but on the other hand 
presently available information on the behavior of the 




(24) 




11 



cross-sectional area of an updraft is scanty and somewhat 
conflicting. 

The system of equations (5)» (12), (l4), (l6), (21), 
and (18) or (20) may be solved first for the cloud 
temperature lapse 

^ j 1 r ) - 1 " c. 



(25) 



d i 






( 26 ) 



dw 

ci H 

(27) Xi 



±1 f 

~ x\ \ 



v'e 



U+Jp) 





V 


i-Jl)- 

'VB 


- 

vV 


.) + qx(T- 




ual 


U- 


t- 








H- i^.Sq) 


r (4t- 4e) + 










L 


the other 


parameters 


A 
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\ JX' 


/ VN/ 






X dz 


)di)^ _ . 






-T 

f v'e I 


cTi 




w \ 





_ \ 



(28) 



( 29 a) 



d ? 

di 



^ ^ _±3 

R dr R,T„ 

- _ (■'+'1-'^,) dn 

d 2 Az 



- hi ) 






or 



(29b) ^ = _ _ A 

d z d z 



W 



Equations (21) through (24) are used as needed. Given 
initial values of T, N, w, jj? » ^ ♦ q » foregoing 

system of equations may be integrated by numerical 
methods over the desired range of z. 
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I ( O iP^ 

* ' ' fi 

fU t i •* » . ^ 








Since adequate data concerning the vertical velocity 
of the entrained environmental air is lacking; eq (2?) 
will be further simplified by assuming that = 0 . 
Moreover since typical values of J? are on the order of 
0OO5 we may place l/fl-i-A) = 1 , The error Introduced 

by the latter approximation will be on the order of ^ 
of one percent. 

If eq (27) is modified as indicated above^ it becomes 



a form which differs from the corresponding equation 
in the Haltiner model only in the last term, which has 
the form 



Similarly in the Haltiner equations for dX/d^r and 



in the modified model. Equations (28) and (29 a) are 
exactly the same as their counterparts in the Haltiner 
model. The only other difference in the two models is 
the appearance of the additional eq (29b). 

4, Computations and Results. 

The foregoing system of equations was Integrated 
by numerical methods on an electronic digital computer, 
the National Cash Register 102A. With a given set of 
initial values, properly scaled inputs to the integration 
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routine were computed for the model given above. These 
values were then integrated over an interval of 10 meters 
yielding a new set of "initial" values for the level z + 

10 meters which served as an input for the computation 
of a new set of Inputs for the integration routine. 
Computation proceeded in this stepwise fashion until the 
top of the cloud was reached. 

Computations were made for four sets of initial 
conditions. In each case the Initial values corresponded 
to one of the sets of initial values used by Haltiner^H^ . 
In Case (a) and Case (b) the same set of initial values 
were used, but in Case (a) only the effect of detrainment 
was considered. The effect of the drag of the liquid 
water was eliminated from Case (a) by an appropriate 
modification to the complete program. In Cases (b) , 

(c), and (d), both the effects of detrainment and the 
drag of the liquid water are Included. 

The initial values used in the computation are as 
follows; 

Case (a) 

Tq = 20 C, Teo ■ 19 °C, Wq = 1 meter/second 
~ ,001/second, ^e = ,7 **0/100 meters 
r " 80 % 

Case (b) 

Same as Case (a) except that both detrainment and 
water drag are Included. 
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Case (c) 

Same as Case (b) except that .6° C/100 meters 

Case (d) 

Same as Case (b) except that r = 90^ 

The results of the computations are presented in 
figures 1 through 6 and table 1. 

Figures 1 and 2 present the results of Cases (a) 
and (b), in which the parameters W, ^ , M/Mq ^ ^ 

are plotted as functions of height. The parameter M/Wq 
is the cumulative change in mass experienced by a unit 
of air due to entrainment. A value of two, for example, 
indicates that the mass has doubled. The parameter AT 
is defined as T - T^ . 

Figures 3 , 5 , and 6 present a comparison of the 

results obtained by Haltlner f or the set of values 

listed as Case (a) with the results of Cases (a) and (b). 
Each parameter Is plotted separately as a function of 
height. The values computed by Haltiner and the values 
from Case (a) should be the same at and below the level 
where detrainment begins. The systematic difference in 
two values below this point is the result of minor 
differences in the computational procedures (primarily 
in the computation of the environmental parameters) 
followed in this investigation and those employed by 
Haltiner. This systematic difference is not significant 
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in itself, but it does illustrate how small differences 
at one level propagate upward and tend to accumulate. 

The effects of detrainment and the drag of the 
liquid water are best illustrated in figures 4 and 5* 

It is apparent that the effect of detrainment alone is 
to accelerate the vertical current. The cloud height 
is greater in Case (a) than it is in either Case (b) or 
the corresponding case computed by Haltlner\_H] , It is 
also apparent from the figures that in Case (a) the 
magnitude of the maximum values of w and M/M^and the 
elevation at which these maxima occur are greater than 
either the values obtained by Haltiner or the values 
of Case (b)o These results, and the results of Cases 
(c) and (d) are summarized in table one, 

5o Conclusions 

The number of cases investigated was not exhaustive, 
but was sufficient to show that the modified model 
yields results that conform quite closely to the results 
of Haltiner[4-1 • In particular, the cloud heights were 
substantially the same but the vertical velocity and the 
amount of entrainment are reduced. These results are 
the natural consequence of the combination of detrainment 
and the drag of the liquid water. Detrainment on the 
one hand provides an acceleration, while the drag of the 
liquid water acts as a retarding force. These forces 
tend to balance one another, but the drag of the liquid 
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water somewhat overbalances the acceleration due to 
detrainment . 

Haltiner investigated an actual cloud for which 

observational data was available, and found that his 
results were in good agreement with the observed values. 
He did find, however, that his values of vertical 
ve loci ty^ although of the right order of magnitude, were 
somewhat in excess of the observed values. To the 
extent that the modified model reduces vertical velocity, 
the modifications appear to be a correction in the right 
direction. 
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Figure 5« Vertical velocity . w, as a fimction of height for 
Case (a), Case (h) and after Ilaltiner. Tide marks on 
curves indicate level at wJiich detrainment stai'ted. 
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